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Abstract 

We derive a mathematical model of progressive collapse and examine role of compaction. Con- 
trary to a previous result by Bazant and Verdure, J. Engr. Mech. ASCE 133 (2006) 308, we find 
that compaction slows down the avalanche by effectively increasing the resistive force. We compare 
currently available estimates of the resistive force, that of Bazant and Verdure (2006) corrected 
for compaction for World Trade Center (WTC) 2, and of Beck, www.arxiv.org jphysics/0609105 , 
for WTC 1 and 2. We concentrate on a damage wave propagating through the building before the 
avalanche that figures in both models: an implicit heat wave that reduces the resistive force of the 
building by 60% in Bazant and Verdure (2006), or a wave of massive destruction that reduces the 
resistive force by 75% in Beck (2006). We show that the avalanche cannot supply the energy to 
the heat wave as this increases the resistive force by two orders of magnitude. We thus reaffirm 
the conclusion of Beck (2006) that the avalanche is initiated in the wake of the damage wave. 
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Bazant and Verdure [l| proposed the following mathematical model to describe the pro- 
gressive collapse in a tall building of a homogeneous longitudinal density po, 

Po - ■^ix2X2) = Rix2) + pogx2, (1) 

where X2 is the position of the avalanche front, po = M/H with M the total mass and H 
the total height of the building and g is the gravity, while R = R{x2) is a local resistive 
force. Here, the dot above the quantity indicates its differentiation with respect to time. 
The compaction ratio n is defined as 

(2) 

where p is a density of the "compacted" section of the building, cf. Fig.[Tl 

We reexamine the steps that lead to Eq. ([1]) from the point of view of classical mechanics. 
As can be seen from Fig.[T] we can choose between two generalized coordinates. The first 
choice is Xi = Xi{t) which for t > well describes the motion of the avalanche. The second 
choice is X2 = X2{t), which represents the position of the avalanche front - an idealized 
point-like boundary between the stationary and the moving part of the building, at which 
the compaction takes place. Motion of the avalanche front is more complex than the motion 
of the avalanche as it combines the motion of the avalanche with its spatial growth due 
to non-zero compaction ratio. While care must be exercised when deriving an equation 
of motion for each of them, the final result may not depend on the choice of generalized 
coordinate. 

For simplicity, we assume that the total energy and the total mass in the system building- 
avalanche is conserved. We note that in that case we obtain the fastest avalanche, as then 
there are no conversion losses of the potential energy of the building into the kinetic and 
then into the crushing energy of the avalanche. Also, the conservation of energy allows us 
to use Lagrangian formalism to derive the equation of motion. 

First, we state the two constraints of descent, 

Xi — Xq = h, (3a) 

PqH = pq {xi - xo) + p (x2 - xi) + Po {H - X2) . (3b) 

Here, the top section stretches from xq = Xo{t) down to Xi = Xi(t), while the compacted 
building occupies from xi down to X2 = X2{t). The point X2 is the avalanche front. With 
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Eq. (l3a|) we state that the length of the top section h does not change in descent, while with 
Eq. fl3bp we express the conservation of the mass of the building. Differentiation of Eq. ([3]) 
with respect to time yields dynamical constraints, and Xi = (1 — k) ±2- 

Second, we find kinetic, potential and latent energy necessary for the Lagrangian formula- 
tion. The kinetic energy K is given hy K = ^ j dx p{x) f ^(x), where the velocity distribution 
in the avalanche is v{x) = xi for x G [xo,X2)- The kinetic energy of the avalanche is thus 

K = ^PoX2xl = ^po{l - K,fx2xl. (4) 

The potential energy U of the whole building is U = — f dx p{x) x g , yielding 

U = -^gpo {H' + {l-K)ixl-h')). (5) 

The latent energy L = L{x2) = — dx' R{x'), produces the resistive force of the building, 
-^(2^2) = -dL{x2) / dx2. 

Given a Lagrangian C = K — U — L, which is a function of a generalized coordinate x 
and its generalized velocity x, the equation of motion follows from ^ dC/dx = dC/dx. We 
recall that we have two choices for the generalized coordinate: x = Xi for the motion of 
avalanche, 01 x = X2 for the avalanche front. If X2 is chosen as a generalized coordinate, we 
simply obtain, 

po (1 - fi:)^^ {X2 X2) = R{x2) + pogil- k)x2 + ^ (1 - x\. (6) 
With x\ as a generalized coordinate we note that d/dxi = (1 — K)^^d / 8x2, yielding 

As expected, the two equations of motion are identical. In Eqs. ([6]) and we introduced 
an additional parameter t] which may take values 1, if the total energy is conserved, or 0, 
if this is not the case. A distinction between the two cases is discussed in Comparison 
between Eq. ((Tj) and Eq. ([1]) of Bazant and Verdure's shows that due to non-zero compaction 
the avalanche propagates through the building faster than it travels, which leads to an 
amplification of the building's resistive force by the factor (1 — > 1. In other words, 
the avalanche front in Bazant and Verdure's model, Eq. ([T]), is faster than the one proposed 
in Eq. ([7]). Eq. ([7]) as a proposed correction to Eq. ([1]) is a major result of this technical note. 
We next discuss the estimates for the resistive force R that can be found in the literature. 
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Bazant and Verdure [l| analyzed collapse of World Trade Center 2 in terms of a model ([T]). 
They assumed that i? is a constant throughout the building's primary and secondary zone, 
where R = AL/AH. For a crushing energy they made an educated guess, AL = 2.4 GNm, 
with the floor height being AH = 3.7 m. Considering the total mass of the building to be 
M = 3.2 ■ 10^ kg, this yields for the resistive force R/{M g) = r = 0.2, as their initial esti- 
mate. They noted that in order for the avalanche to reach the ground level in 2T ~ 10.8 s 6| 
they had to use R/2 instead of R. The 50% reduction, they argued, came from heat (p. 15, 
top paragraph of the on-line edition of jl|). That is, an assumption built in their model is 
that the avalanche pushes a heat wave in front of itself which reduces the strength of the 
building by 50%. Using their value for compaction k, ~ 0.2 and the corrected equation of 
motion proposed here, requires R to be reduced by 60% instead. 

Here the following comment is in place. If the avalanche were supplying energy to the heat 
wave then R in Eq. ([1]) splits in two components. First, original R decreases, say, by 60% 
as discussed earlier. However, an additional resistive force {R)heat appears which describes 
the rate of transfer of energy from the avalanche to the heat wave per unit length. A 
simple estimate shows [7] that {R)heat/ {M g) ~ 20 ^ i?. It is thus obvious that the 60%- 
strength-reducing heat wave could not have been created or maintained by the avalanche. 
On the contrary, if the avalanche were in fact heating the steel, then this acted as a resis- 
tive force comparable to R. E.g., the avalanche warming the steel by AT = 10 K yields 
{R)heat/{M g) = r 0.1 which is comparable to the observed resistive force. We do men- 
tion, however, that the authoritative document on collapse of World Trade Center 1 and 2, 
the NIST report |3i] exphcitely states that no elevated temperatures were observed in the 
secondary zone. 

Rather than guessing, we proposed a procedure for estimating -R {4], where one first finds 
the ultimate yield strength Y of the vertical columns using their specifications, following 
which the resistive force R is estimated from a simple linear model, R = e-Y , with e = 0.25 
being the ultimate yield strain of structural steel under compression. Applying this to 
WTC 1 and 2 led to an initial estimate of R{x2)/{M g) = r + s ■ {X2/H), with r ^ 0.2 
and s ^ 0.7. We analyzed the collapse of World Trade Centers 1 and 2 where we divided 
the building into the primary and the secondary zone where we allowed Ri in the primary 
zone to be considerably smaller than Ru in the secondary zone, Ri = 1/4 ■ Rjj for WTC 1 
and Rj = 3/8 ■ Ru for WTC 2, while we neglected the compaction altogether. We did this 
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subdivision after the statement from the NIST report that the damage to the buildings was 
concentrated in their primary zones while leaving the secondary zones intact. We found that 
to reach the collapse time 2T the initial estimate of R had to be reduced by 75%. In our 
report we dubbed the 75%-strength-reducing wave that preceeded the avalanche the wave of 
massive destruction (WMD). As the avalanche was not supplying the energy to the WMD, 
and the WMD propagated before the avalanche, we concluded that the WMD caused the 
avalanche. 

Currently, we can only speculate about the source of the 60- 75%-strength-reducing wave 
and its coupling to the avalanche. However, a piece of information that would provide an 
important insight is the descent curve, which describes position as a function of time of some 
visible part of the building, say its top, Xq = Xo{t). Once the descent curve is known it is the 
acceleration, = xo{t), that can be directly connected to R through a mathematical model. 
In fact, in 5| we examine the descent curve available for WTC 7 in terms of a corrected 
model of "crush- up" mode of progressive collapse and identify the phases of descent. 

Lastly, as the collapses of World Trade Centers resemble controlled demolition it would 
be instructive to apply the methods discussed in this and other articles to other buildings 
that are known to have been destroyed in controlled fashion. 
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at which steel strength is reduced by 50%, and 0.5 < cr < 0.9 for the ratio of the mass of steel 
to the total mass of the building. This yields {R)heat/{M 9) = crC AT/{H g) = 17.8, which is 
by two orders of magnitude greater than r ^ 0.2, the apparent resistive force of the buildings. 
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FIG. 1: Propagation of an avalanche in a tall building of uniform density pQ. Extraneous factors 
cause initial weakening of the load bearing structure in the primary zone, leaving the building below 
(secondary zone) intact. The avalanche forms at the top of the primary zone, which propagates 
and compacts the building underneath from po to p, with k = po/ p <^ 1 being a compaction ratio. 
The idealized scenario allows us to identify the following points: xq, the top of the building; xi, 
beginning of the compacted section of the avalanche; and X2, the location of the avalanche front. 
Here it is implied that it is easier for the avalanche to drop ("crush down") then to stop and 
compact the top section instead ("crush up"). 
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